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Detecting metal-rich intermediate-age Globular Clusters 
in NGC 4570 using K-band photometry 
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Abstract Globular Cluster Systems (GCSs) of most 
early-type galaxies feature two peaks in their optical 
colour distributions. Blue-peak GCs are believed to be 
old and metal-poor, whereas the ages, metallicities, and 
the origin of the red-peak GCs are still being debated. 
We obtained deep K-band photometry and combined it 
with HST observations in g and z to yield a full SED 
from optical to near-infrared. This now allows us to 
break the age-metallicity degeneracy. 

We used our evolutionary synthesis models GALEV 
for star clusters to compute a large grid of models with 
different metallicities and a wide range of ages. Com- 
paring these models to our observations revealed a large 
population of intermediate-age (1-3 Gyr) and metal- 
rich (« solar metallicity) globular clusters, that will 
give us further insights into the formation history of 
this galaxy. 



1 Introduction 

Globular Cluster Systems (GCSs) are now recog- 
nized as powerfu l tracers of their p arent galaxy's for- 
mation history dWest et al.l 12004 lFritze-v. A.I 12004 ; 
Brodie and Straderll2006ir From their age and metallic- 
ity distributions one can reconstruct the parent galaxy's 
(violent) star formation and chemical enrichment his- 
tory all the way from the very onset of star formation 
in the early universe to the present. 

Most early-type galaxie s show a bimodal color distri- 
bution for their GCSs (e.g. Gebhardt and Kissler-Patigl 
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19991: iKundu and WhitmOTell200ll |?l; IPeng et alJEooi ): 
A universal blue peak and a red peak with colors and 
height relative to the blue peak varying from galaxy to 
galaxy. The blue peak GCs are believed to be old and 
metal-poor, the origin of the red peak is still unclear. 

1.1 Our approach to lift the age-metallicity 
degeneracy 

Optical data alone do not allow to disentangle ages 
and metallicities: Colour-to- metallicity transforma- 
tions have to assume an age, while colour-to-age trans- 
formations are only valid for one metallicity. This de- 
generacy, however, can be broken by including near- 
infrared data that are more sensitive to changes of 
m etallicity rathe r than age. 



Anders et al.l (|2004l ) used extensive artificial star 



cluster tests and showed that observations in 



3 passbands for GCSs in dust-free E/SOs or 

4 passbands for (young) star clusters in dusty envi- 
ronments, 

spanning as wide as possible a wavelength-basis (U 
through K) and 

including at least one NIR-band (e.g. K) with 
accuracies < 0.05 mag in the optical and < 0.1 mag 
in the NIR 



allow to disentangle ages and metallicities and deter- 
mine individual GC metallicities to < ±0.2 dex, 
and ages to < ±0.3 dex, i.e. they allow to distinguish 
< 7 Gyr old GCs from those > 13 Gyr old. 

Similar studies also using NIR-data to determine 



been done for a few gal 


ixies until now (Puzia et al. 


20021 Kissler-Patig et al. 


2002; 


Hemoel et al. 2003: 


Larsen et al.l 2005: Hcmpel et al. 


2007h. More than 
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half of these are found to host a population of GCs 
that is younger and/or more metal-rich than the old 
and metal-poor GC population in the Milky Way. 
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1.2 1.4 
(F475W - F850LP) 



Fig. 1 Color-color diagram showing evolutionary 
tracks of globular clusters of different metallicity in the 
(F850LP - K s ) - (F475W - F850LP) plane, computed 
from our GALEV evolutionary synthesis models and using 
AB magnitudes; Tick-marks are shown every 2 Gyrs. The 
vertical lines indicate the peak-colo urs of the bimo dal 
optical colour distribution (taken from lPeng et al.ll2006l ). 



We note that many of these previous studies relied on 
cumulative distribution. The most important change 
compared to our analysis is that we derive the physical 
parameters for each individual cluster. That also en- 
ables us to study correlations of these parameters, e.g. 
with their spatial distributions. 



2 Models 

We used our GALEV evolutionary synthesis models for 
star c lusters ( Schulz et al.ll2002 : Anders and Fritze-v. A.^ 
20031 ) to compute a large grid of models for five differ- 
ent mctallicities —1.7 < [Fe/H] < +0.4 and ages be- 
tween 4 Myr and 16 Gyr with time-steps of 4 Myr. Since 
early-type galaxies do not contain significant amounts 
of dust we did not include extinctions E(B — V) > 
into our grid. We therefore only need three filters (HST 
F475W,F850LP and SOFI Ks) to determine all relevant 
parameters (age, metallicity, mass) for each cluster. 

Note that we do not depend on color-transformation 
from the HST to Standard Johnson filters. Our models 
first compute spectra as function of time, that later are 
convolved with the corresponding filter curves to yield 
final magnitudes. Figure [1] shows the resulting color- 
color diagram for the three filters discussed here. 

As expected, there is a degeneracy between ages and 
metallicities, but note how nicely the addition of the 
near-infrared splits up different metallicities and allows 
to break this degeneracy. 



3 Data 

3.1 Near-infrared observations 

We used the SOFI (Son of ISAAC) near-infrared im- 
ager on the ESO-NTT during two nights in May 2007 
(ESO-program id 079.B-0511). SOFI contains a Hawaii 
HgCdTe 1024x1024 chip with a resolution of 0.288 
arcsec per pixel, resulting in a field-of-view (FoV) of 
4.92 x 4.92 arcmin. This large FoV allowed us to sig- 
nificantly reduce the overhead for sky-exposures using 
one half of the detector for the galaxy, the other half 
for sky, and changing positions every minute. 

3.2 K-band data reduction 

The data reduction was done using ESO-MIDAS and 
largely followed the r ecipes given in the SOFI User's 
manual <|Sterzikll2007h . starting with the inter-quadrant 
row-crosstalk. To flat-field the data we used a com- 
bination of dome-flats, illumination correction surfaces 
and a refined master-flat obtained from normalized sky- 
frames of both nights. All frames were sky-subtracted 
using the average of six frames taken nearest in time 
that have been scaled to the sky-value of the object 
exposure. We aligned individual images by matching 
the positions of several stars and compact background 
galaxies and averaged them to give the final image. 

3.3 HST data 

We used archival data obtained from the Hubble Space 
Telescope. Both datasets (J8FS18011 and J8FS18021) 
were observed as part of The ACS Virg o Cluster Sur- 
vey ( Cote et "all 2004 ; Jordan et al.ll2004 ) and automat- 
ically reduced and calibrated by the On-the-fly Repro- 
cessing (OTFR) pipeline at STScI. After retrieval we 
checked the alignment of both frames relative to each 
other, again using a set of stars within the ACS FoV. 

3.4 Cluster selection and photometry 



Clust er selection was done using SExtractor ( Bertin and Arnouts^ 
1996h requiring at least 4 pixels with intensities above 
the 3<r-over-background threshold. A cross-correlation 
of the catalogs of both HST filters to remove remain- 
ing spurious detections resulted in a "optical" cata- 
log of 330 sources. For all these cluster candidates 
we derived intrins ic radii using the ISHAPE package 
within BALOAB (iLarsenl 119991 ), assuming a circular- 
symmetric King profile with concentration c = 30 for 
the GCs. 

To remove background galaxies and stars, we re- 
quired all valid cluster candidates to have radii in 
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Fig. 2 Distribution of ages (upper panel) and metallicities 
(lower panel) for all clusters 



the range 0.2px < r < 5px (equivalent to physical 
sizes of 0.8 . . . 20 p c at the assumed distance of 17Mpc 
( Tonrv et aT1l200lh ), leaving w 280 candidates. For all 
of them we performed aperture photometry with aper- 
ture sizes of 10 pixels for the HST filters and 7 pixels 
(» 2 arcsec) in Ks. This allowed us to obtain (g, z, Ks) 
photometry for rs 150 candidates; most of the remain- 
ing cluster candidates for which we could not derive 
magnitudes in all three filters were not included within 
the SOFI FoV. 



4 Results 



We derived physic al parameters for a ll GC candidates 



using AnalySED (jAnders et all [2004,?). It 



compares 

the observed spectral energy distributions (SEDs) with 
all model SEDs and automatically derives probabilities 
for each SED on the basis of a x 2 -algorithm. From 
these probabilities it finds the best-fitting template and 
its physical parameters age and metallicity. The results 
are shown in Figure [5J 

4.1 Old and intermediate ages 

The age distribution is dominated by two distinct popu- 
lations: An old population with ages older than 10 Gyr 
that has formed during galaxy formation in the early 
universe. This population is universal in all galaxies 
and can be studied in great detail within our Milky 



Way. The second population has intermediate ages of 
1-3 Gyr and therefore must have formed later during 
a violent episode of star-formation; the most plausible 
explanation for such an event being an intense starburst 
as e.g. accompanying the merger of two gas-rich spirals 
or the accretion of a gas-rich companion, resulting in a 
phase of massive star cluster formation. Besides these 
two dominating populations there is a number of clus- 
ters that do not belong to either of these two. These can 
be explained by gas left over from the original galaxies 
that was ejected into tidal tails and later rained down 
onto the merger remnant. 

A merger or accretion event is supported by the 
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1998 
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(1998) estimated an age of < 2 Gyr for the central struc- 
ture, in excellent agreement with our ages for the inter- 
mediate population, suggesting that both have formed 
from the same event. 

4.2 Solar metallicities 

The lower panel of Fig. [2] shows the metallicity dis- 
tribution of our cluster sample. It features an eye- 
catching prominent peak at solar metallicity. 26 out 
of 36 of these high-metallicity clusters belong to the 
intermediate-age population, however, 10 out of the 36 
solar metallicity GCs have old ages. Their origin is not 
yet clear. 

We caution the reader that the number of old and 
metal-poor globular clusters is underestimated with our 
analysis. Since we only include globular clusters with a 
secure Ks-band detection, we are biased toward higher- 
metallicity clusters, because those have significantly 
redder optical-near-infrared colours. 

A more detailed analysis of the properties of our clus- 
ter sample, including their masses will be published in 
Kotulla et al., in prep. 



5 Summary 

We obtained deep K-band photometry of the Virgo 
lenticular NGC4570 and combined it with archival opti- 
cal data from HST. We selected globular cluster candi- 
dates based on the HST data and their intrinsic optical 
sizes. 

GALEV evolutionary synthesis models were used in 
combination with AnalySED to automatically derive 
physical parameters age, metallicity and mass for each 
individual cluster. 

We detect a significant population of intermediate- 
age (1-3 Gyr) and metal-rich ([Fe/H] > —0.4) cluster 
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population that has not been reported for this galaxy 
before. 
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